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DIGEST 


v 

_ / 

A  vorticity  component  and  budget  analysis  in  three- 
dimensions  is  conducted  for  the  subcloud  layer  of  two 
microburst-producing  storms  using  dual-Doppler  derived  winds. 
The  raw  data  were  collected  during  the  Joint  Airport  Weather 
Studies  (JAWS)  project."  This  research  examines  the  vorticity 
associated  with  a  single  microburst-producing  storm  which 
occurred  on  14  July  1982  and  a  multiple  microburst-producing 
storm  which  occurred  on  5  August  1982.  The  calculated  vorti¬ 
city  data  are  examined  over  the  storm  domain  (10  km  x  10  km  for 
the  14  July  case  and  15  km  x  15  km  for  the  5  August  case)  and 
the  microburst  domain  (5  km  x  5  km)  including  two  analysis 
times  for  the  first  case  and  three  for  the  latter  case. 

Results  indicate  that  horizontal  vorticity  centers  are 
coincident  with  strong  horizontal  gradients  of  vertical  velo¬ 
city  and  with  areas  of  strong  vertical  shear.  These  gradients 
are  maximized  along  the  edges  of  downrushing  air  and  along  gust 
fronts.  Regions  of  high  speed  low-level  winds  have  a  core  of 
horizontal  vorticity  above  them. 

A  vorticity  budget  analysis  of  the  advection,  divergence, 
and  tilting  terms,  comprising  the  vorticity  component  equa¬ 
tions,  shows  the  magnitude  of  these  terms  to  be  greater  in  the 
microburst  domain  than  in  the  storm  domain.  This  indicates  the 
strongest  forcing  and  advection  takes  place  within  small 


1 


2 


regions  of  the  storms.  The  diverging  outflow  of  the  microburst 
in  both  cases  weakens  any  existing  positive  vorticity  in  the 
microburst  region,  e.g.,  below  0.75  km. 

The  differences  in  storm  structure  allow  examination  of 
the  vorticity  of  a  simple,  nearly  circular-symmetric  microburst 
and  that  of  a  microburst  within  a  complicated  flow  field.  This 
leads  to  different  vorticity  distributions  and  budgets.  )(< 
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1.  Introduction 


The  first  measurements  of  airflow  associated  with  micro¬ 
bursts  were  made  with  Doppler  radar  during  the  1978  Northern 
Illinois  Meteorological  Research  on  Downbursts  (NIMROD)  Pro¬ 
ject.  These  measurements  can  only  be  treated  as  single  Doppler 
measurements  because  the  large  spacing  between  radars  (60  km 
baseline)  and  the  curvature  of  the  earth  caused  resolution 
problems,  i.e.,  obstruction  of  low-level  wind  observations  by 
Doppler  radars,  past  30  km.  NIMROD  researchers,  however,  docu¬ 
mented  10  microbursts,  5  gust  fronts,  and  two  supercells,  but 
again,  the  main  difficulties  were  in  the  ground  clutter  and 
curvature  of  the  earth. 

The  Joint  Airport  Weather  Studies  (JAWS)  program  conducted 
research  through  the  summer  of  1982  to  determine  the  fine-scale 
structure  of  thunderstorm  kinematics  in  the  vicinity  of 
Denver's  Stapleton  International  Airport  (Fujita,  1985).  The 
network  included  a  triple  Doppler  radar  (CP-2,  CP-3,  and  CP-4) 
triangular  network  (maximum  baseline  of  28  km),  27  Portable 
Automated  Mesonet  (PAM)  stations,  1  Low-Level  Wind  Shear  Alert¬ 
ing  System  (LLWSAS),  2  CC^  lidar  radars,  and  5  research  air¬ 
craft  (Figure  1.1). 

Fujita’s  (1985)  definition  of  a  downburst  as  a  strong 
downdraft  which  induces  an  outburst  of  damaging  winds  on  or 
near  the  ground  expanded  as  a  result  of  NIMROD  to  include  two 
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Fig.  1.1:  Map  illustrating  JAWS  Project  facilities 
situated  in  the  vicinity  of  Denver's  Stapleton 
International  Airport  (after  McCarthy  et  al.,  1983). 
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sub-categories,  the  raacroburst  and  the  microburst.  This  was  to 
account  for  the  varied  length  scales  of  observed  downbursts. 
The  microburst,  then,  became  defined  as  a  strong  downburst  with 
outflow  diameter  of  less  than  or  equal  to  4  km  and  peak  winds 
lasting  2  to  5  minutes. 

Wilson  et  al.  (1984)  modified  Fujita's  definition 
slightly  to  refer  to  Doppler  radar-observed  diverging  outflows 
near  the  surface  associated  with  convective  storms  during  JAWS. 
The  differential  Doppler  velocity  across  the  divergence  center 
must  be  greater  than  or  equal  to  10  m/s  with  initial  distance 
between  approaching  and  receding  maximum  velocity  centers  being 
less  than  or  equal  to  4  km.  Using  this  criterion,  about  70 
microbursts  were  identified. 

These  intense,  divergent  outflows  at  small  scales  are  con¬ 
sidered  to  be  very  dangerous  to  aviation  safety.  In  addition 
to  diverging  outflow,  horizontal  vortices  which  are  common 
features  of  microburst  outflows  (Hjelmfelt,  1987)  have  been 
suggested  as  the  cause  of  1985  Dallas-Fort  Worth  L-1011  jet 
aircraft  crash  (Fujita,  1986). 

Figure  1.2  shows  a  schematic  in  the  vertical  of  an  evolv¬ 
ing  microburst  wind  field  (Wilson  et  al.,  1984).  Notice  at 
time  T-2  minutes  the  downdraft  begins  to  diverge  horizontally 
and  lower  to  below  1  km  with  possible  horizontal  roll  vortices 
developing  on  the  edges  of  the  outflow.  With  time  the  vortex 
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descends  further  and  the  surface  winds  increase.  This  figure 
demonstrates  the  rapid  development  and  evolution  typically 
found  with  such  events. 

Kessinger  et  al.  (1988)  studied  the  sub-cloud  layer  of  a 
multicellular  storm  in  Colorado  using  multiple  Doppler  derived 
winds.  Their  storm  produced  misocyclones  which  are  horizontal 
cyclonic  circulation  centers  with  diameters  of  2-4  km  (Fujita, 
1985),  downbursts,  and  horizontal  vortex  circulations.  Many 
studies  have  examined  the  rotation  (vorticity)  generation  as  it 
pertains  to  the  mesocyclones  that  produce  tornadoes.  Kessinger 
et  al.  (1988)  found  that  the  misocyclone  characteristics, 
differ  from  those  of  mesocyclones  by  having  a  vorticity  maximum 
near  cloud  base  instead  of  at  low  levels  and  that  the  low-level 
positive  vertical  vorticity  is  weakened  by  the  low-level  diver¬ 
gence  associated  with  the  raicroburst.  They  also  point  out  that 
horizontal  vortex  circulations  (rotors)  form  along  the  edge  of 
these  misocyclones  and  storm  downdrafts  and  propagate  away  from 
the  storm.  These  rotors  have  also  been  associated  with  the 
regions  of  maximum  surface  winds. 

The  importance  of  the  rotor  in  causing  damage  as  well  as 
being  an  aviation  hazard  has  recently  come  under  more  study. 
The  rotor  core  from  the  microburst  outflows  has  been  associated 
with  lower  pressure  than  its  surroundings,  which  acts  to 
accelerate  the  surface  winds  (Waranauskas ,  1985)  leading  to  an 
explanation  for  burst  swath  damage  and  why  surface  wind  speeds 
are  greater  than  the  microburst  downflow  wind  speeds. 
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Waranauskas  (1985)  further  suggests  that  the  axis  of  the  rotor 
and  the  microburst  coincide  thereby  linking  the  rotor  as  a 
cause  or  enhancement  of  the  microburst.  Linden  and  Simpson 
(1985)  state  that  the  wind  shear  and  downward  motion  associated 
with  the  back  of  the  vortex  (rotor)  may  be  responsible  for  the 
danger  of  flying  through  a  microburst.  An  additional  mechanism 
for  intensifying  the  leading-edge  vortex  may  be  the  existence 
of  rotation  in  the  descending  air. 

Despite  the  importance  of  the  microburst,  misocyclones , 
rotors,  etc.,  few  studies  have  been  conducted  outside  the  con¬ 
cern  for  the  vertical  component  of  vorticity.  This  research 
will  utilize  the  derived  wind  fields  from  dual-Doppler  radar 
data  by  Lin  et  al.  (1987),  Lin  and  Hughes  (1987),  and  Lin  and 
Coover  (1988)  to  study  the  vorticity  (in  three  dimensions) 
associated  with  microburst-producing  storms. 


2.  Statement  of  Problem 


The  main  objective  of  this  research  is  to  study  the  vorti- 
city  associated  with  raicrobursts.  Since  microbursts  have  rota¬ 
tion  around  horizontal  (e.g.,  rotors)  and  vertical  (e.g.,  miso- 
cyclones)  axes,  a  three-dimensional  analysis  is  required. 

Therefore,  a  three-dimensional  vorticity  budget  analysis 
for  the  subcloud  layer  of  two  microburst-producing  storms  will 
be  examined  using  wind  fields  derived  from  JAWS  dual-Doppler 
radar  data.  The  three-dimensional  vorticity  equation  (Appen¬ 
dix)  in  vector  form  will  be  separated  into  its  three  scalar 
component  equations.  The  specific  terms  being  investigated  are 
the:  1)  horizontal  advection;  2)  vertical  advection;  3)  vorti¬ 
city  divergence;  and  4)  tilting  terms. 

The  first  storm  which  occurred  on  14  July  1982,  produced  a 
single  microburst,  a  rotor,  and  a  gust  front.  Analysis  is  per¬ 
formed  for  1647  and  1649  MDT.  The  storm  was  at  its  quasi¬ 
steady,  mature  stage  at  1647  MDT  and  began  to  decay  by  1649 
MDT.  The  vorticity  vector  remained  quasi-horizontal  through 
both  time  periods.  Budget  information  should  help  determine 
why  no  misocyclone  was  produced  in  this  case. 

The  second  storm,  which  occurred  on  5  August  1982,  pro¬ 
duced  two  microbursts,  a  microburst  misocyclone,  and  a  gust 
front.  Analysis  is  performed  for  three  time  periods,  1845, 
1847,  and  1850  MDT.  This  storm  will  allow  the  study  of  raisocy- 
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clone  evolution. 

Finally,  a  comparison  can  be  made  as  to  the  mechanisms 
responsible  for  one  storm's  vorticity  field  to  remain  quasi¬ 
horizontal  while  another  developed  rotation  around  the  vertical 
axis.  The  ultimate  goal  is  to  understand  the  importance  of  the 
microburst  and  its  associated  structure  to  the  generation  of 
vorticity  through  interpretation  of  the  various  terms  of  those 
component  equations  comprising  the  three-dimensional  vorticity 
equation. 


3.  Methodology 

3.1  Storm  Specifics 

The  three  wind  components  used  to  conduct  this  research 
were  derived  from  dual-Doppler  radar  data  (Lin  et  al.,  1987; 
Lin  and  Hughes,  1987;  Lin  and  Coover,  1988).  The  data  were 
collected  in  Colorado  during  JAWS  in  the  summer  of  1982.  The 
first  case  is  the  14  July  1982  case  (or  the  simple  case)  from 
which  a  single  microburst  was  produced  along  with  a  gust  front 
and  horizontal  vortex  (rotor)  circulation  (Hughes,  1986; 
Coover,  1988).  This  case  may  be  categorized  as  a  quasi¬ 
stationary,  single  microburst  event.  As  such,  the  specific 
structure  of  microburst  vorticity  can  be  isolated  without  con¬ 
tamination  from  interactions  with  other  raicroburst  outflows. 
This  storm  consists  of  a  warm  core  downburst  of  air  below  0.75 
km  within  an  area  of  moderate  reflectivity  values.  The  outflow 
pattern  is  nearly  circular  symmetric.  There  are  two  analysis 
time  periods,  1647  and  1649  MDT.  At  1647  MDT,  both  radars  were 
able  to  scan  the  entire  storm  volume  (0.25  to  8.5  km),  but  only 
data  at  and  below  2  km  will  be  discussed  since  the  1649 
analysis  period  has  no  data  above  this  level. 

For  the  14  July  1982  case,  the  analysis  grid  is  a  10  km  x 
10  km  grid  centered  on  the  microburst  with  0.5  km  grid  spacing 
in  the  horizontal  and  0.25  km  grid  spacing  in  the  vertical 
below  1  km  and  0.5  km  grid  spacing  above  this  level.  There- 
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fore,  six  levels  of  the  storm  are  considered  at  each  analysis 
time.  This  is  the  storm  domain.  A  sub-domain  or  microburst 
domain  has  been  defined  as  a  5  km  x  5  km  horizontal  grid  cen¬ 
tered  on  the  microburst  with  the  same  grid  spacings  as  in  the 
storm  domain. 


The  second  case  is  the  5  August  1982  case  or  complex  case 
which  produced  two  microbursts,  a  microburst  misocyclone,  a 
gust  front  mesocyclone-like  circulation,  and  gust  front 
(Coover,  1988).  In  this  case,  the  scan  strategy  limited  the 
top  of  the  scanned  volume  to  the  1.25  km  level.  The  analysis 
grid  is  a  15  km  x  15  km  grid  with  a  0.5  km  grid  spacing  hor¬ 
izontally  and  0.25  km  vertical  spacing  beginning  as  with  the  14 
July  case  at  the  0.25  km  level,  but  for  only  five  levels.  The 
microburst  domain  is  a  5  km  x  5  km  analysis  grid  as  defined  for 
the  14  July  case.  Three  analysis  times  are  available  for  the  5 
August  case;  1845  ,  1847,  and  1850  MDT. 


3.2  Vorticlty  Analysis 


Values  for  the  specific  vorticity  components,  £,  n,  and  C» 
were  obtained  by  calculating  the  shear  at  each  grid  point  and 
subtracting  the  respective  values  based  on  the  following  defin¬ 
itions: 


£  ■  9w/9y  -  9v/9z  ;  r|  -  9u/9z  -  9w/9x; 


(1) 


C  -  9v/9x  -  9u/9y . 

£  represents  rotation  in  the  y-z  plane  (about  the  x-axis), 
n  represents  rotation  in  the  x-z  plane  (about  the  y-axis),  and 
C,  most  commonly  used  to  describe  atmospheric  circulations, 
represents  rotation  in  the  x-y  plane  or  horizontal  plane  (about 
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the  z-axis).  A  positive  value  for  each  term  is  defined  as 
cyclonic  rotation  about  each  axis  in  the  positive  axis  direc¬ 
tion  using  a  right-hand  Cartesian  coordinate  system. 

In  a  purely  qualitative  sense,  one  would  expect  signifi¬ 
cant  values  of  £  and  n  in  regions  of  large  horizontal  gradients 
of  vertical  velocity  or  in  regions  of  strong  vertical  shear  of 
the  horizontal  wind,  such  as  those  associated  with  rapidly  mov¬ 
ing  upward  or  downward  flows  in  thunderstorms  or  in  the  lower 
layers  of  a  microburst.  Likewise,  the  C  component  will  be 
small  in  divergent  flow  (irrotational )  and  large  in  regions  of 
strong  horizontal  shear. 

The  vorticity  equation  in  vector  form  is  obtained  by 
taking  the  curl  of  the  equation  of  motion  in  vector  form.  Upon 
recognizing  the  small  effect  the  rotation  of  the  earth  has  on 
the  absolute  vorticity,  20,  at  the  convective  scale,  one 
obtains  an  expression  for  the  relative  (relative  to  the  earth) 
vorticity  vector,  q: 

dq/dt  -  3q/dt  +  V«7q  +  w3q/3z 

^  #  «  »  *  jk  ^  (2) 

-  -q(7h-  Vh+3w/3z)  +  q.VV  -  7  X  (cr7P)  +  7  X  F 
where  q  -  (It  +  Jn  + 

Once  terms  are  expanded,  grouped,  and  the  correction  to 
the  earth  relative  velocity  is  made  (due  to  the  storm  motion), 
one  obtains  the  vorticity  component  equations  in  moving  coordi¬ 
nates  (following  vector  dot  product  operations): 
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3£/3t  -  ~Vh*Vh£  -w3f;/3z  -£(V  •  Vh+3w/3z)  +(^3u/3x+n9u/3y+C3u/3z) 

-(3cr/3y3p/3z-3a/3z3p/3y)  +(3f  /3y-3f  /3z) 

z  y 

^(A)  (B)  ^  £C)  (D) 

3r)/3t  «  -Vh*  Vhn  -w3n/3z  -n(V^*  V^+3w/3z)  +(£3v/3x+n3v/3y+C3v/3z) 

-(3a/3z3p/3x-3cr/3x3p/3z)  +(3f  /3z-3f  /3x)  (3) 

(E)  X  (F)  2 

3C/3t  «  -Vh-  7hC  -w3C/3z  -C(?h*  V^+dw/dz)  +(£3w/3x+n3w/3y+C3w/3z) 

-(3a/3x3p/3y-3a/3y3p/3x)  +(3f  /3x-3f  /3y) 

y  x 

The  result  is  six  terms  per  component  equation  including: 
A)  horizontal  advection;  B)  vertical  advection;  C)  divergence; 
D)  tilting;  E)  solenoidal  term;  and  F)  friction  term.  The 
friction  term  is  generally  one  order  of  magnitude  smaller  than 
the  others  and  will  not  be  examined  here.  Likewise,  the 
solenoid  term  requires  retrieved  information  derived  from  the 
recovered  winds  and  is  subject  to  larger  uncertainty,  hence,  it 
will  not  be  treated  here.  The  other  terms  can  be  determined 
from  the  three  wind  components  and  their  spatial  derivatives. 

A  fourth-order  finite  differencing  scheme  will  be  applied 
in  the  horizontal  and  a  second-order  finite  differencing  scheme 
will  be  used  in  the  vertical  to  obtain  an  estimate  for  each  of 
the  remaining  terms.  The  magnitude  of  these  terms  can  then  be 
determined  by  taking  their  respective  area  means  and  standard 
deviations  over  the  domains  of  interest.  The  vorticity  budget 
in  three  directions  can  be  studied  using  these  results. 
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Therefore,  the  three  component  equations  of  the  3-dimensional 
vorticity  equation  used  in  this  research  are  derived  in  the 
Appendix  but  listed  here  for  convenience: 

3£/3t  -  -V^«  Vh£  -w3£/3z  -£(Vh‘  V^+3w/3z)  +  (£3u/3x+r|3u/3y+C3u/3z) 

_HAD  VAD  ^  DIV  TILT  (4) 

3n/3t  »  _Vh*  ^h11  -w^r'^®z  “h(Vj^  Vh+3w/3z)  +(£3v/3x+n3v/3y+C3v/3z) 

3C/3t  «  -V^»  V^C  -w3C/3z  “C(^h*  V^+Sw/dz)  +(£3w/3x+n3w/3y+C3w/3z) 

HAD  is  the  horizontal  advection  of  the  vorticity  component 
■in  question,  while  VAD  is  simply  the  vertical  advection  of  that 
quantity.  For  a  given  flow  in  a  quasi-steady  state  flow  field, 
these  terms  will  be  largest  in  the  regions  of  tightest  vorti¬ 
city  gradients.  DIV  represents  the  generation  (positive)  of 
the  respective  vorticity  component  through  convergence  and  the 
last  term,  TILT  represents  the  generation  of  vorticity  by  tilt¬ 
ing  the  vorticity  vector  into  the  respective  positive  axial 
direction.  The  solenoidal  term  and  friction  terms  are  not  con¬ 
sidered  here  since  Lin  (1988)  showed  that  friction  is  one  order 
of  magnitude  smaller  than  other  terms  and  the  solenoid  values 
are  subject  to  greater  uncertainties  having  been  calculated 


from  a  derived  wind  field. 


4.  Results 


4.1  Simple  Case:  14  July  1982,  1647  MPT 

Figure  4.1  shows  the  contoured  vorticity  analysis  field 
for  each  vorticity  component  on  the  horizontal  velocity  field 
at  the  0.5  km  level.  The  gust  front  is  denoted  by  the  dashed 
line  to  the  northwest  of  the  center  of  the  microburst,  M.  The 
microburst  domain  is  represented  by  the  small  box  centered  on 
the  microburst.  £  has  positive  values  north  of  the  microburst 
and  negative  to  the  south. 

The  tightest  horizontal  gradient  of  £  appears  across  the  gust 
front  to  the  northwest,  in  evidence  of  the  switch  in  sign  of  £ 
across  the  front  (Fig.  4.1a). 

A  south  to  north  (S-N  on  Figure  4.1)  crossection  of  f;  is 

shown  in  Figure  4.2a.  The  gust  front  is  evident  (dashed  line) 

from  the  velocity  field  north  of  the  raicroburst.  A  rotor  or 

horizontal  vortex  circulation  is  very  apparent  in  the  wind 

field  centered  at  0.75  km  on  the  microburst  side  of  the  gust 

front.  The  £  field  corresponds  to  a  maximum  at  this  location 
_3 

of  18  x  10  per  second.  On  the  south  side  of  the  raicro¬ 
burst,  the  rotor  is  non-existent  although  the  flow  field  does 
show  a  circulation  but  it  is  not  closed.  Values  of  £  are 
strongly  negative  here  but  not  as  significant  as  in  the  rotor 
to  the  north.  The  southern  region  being  on  the  downwind  side 
of  the  storm  may  be  the  reason  the  rotation  is  not  as  evident. 
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Figure  4.1:  Plan  view  of  horizontal  velocity  and  contours  of 
(a).  £,  (b).  n,  and  (c).  C  for  14  July  1982,  1647  MDT  at  0.5  km 
level.  Microburst,  M,  microburst  domain  (small  box),  and  south 
to  north  cross  section  (S-N)  are  also  depicted. 
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That  is,  the  outflow  to  the  north  is  going  against  the 

environmental  flow  while,  to  the  south  the  outflow  is  in  the 

same  direction  as  the  environmental  flow;  therefore,  the  rotor 

does  not  form.  The  £  field  is  closed  in  the  southern  region, 

-3 

however,  with  a  value  greater  than  10  x  10  per  second,  indi¬ 
cating  a  vortex  of  modest  strength  exists  there. 

Values  of  £  are  small  at  the  microburst  center  up  to  at 
least  0.75  km,  these  are  small  compared  to  the  relatively 
higher  values  along  the  periphery  of  the  microburst.  The  flow 
in  this  region  and  in  this  plane  may  then  be  predominantly  non- 
divergent  or  rotational,  whereas  on  the  south  side  of  the 
microburst  the  flow  perhaps  is  composed  of  a  stronger  divergent 
portion  as  well  as  a  rotational  portion.  It  follows  then  from 
this  discussion  that  within  the  microburst  center  the  flow  is 
largely  irrotational  or  divergent  with  respect  to  the  y-z 
plane.  The  divergent  idea  follows  straightforwardly  while  the 
irrotational  aspect  may  not  be  as  obvious.  A  cross-section  of 
vertical  velocity  (Figure  4.3)  shows  the  larger  values  of  posi¬ 
tive  and  negative  £  occur  coincident  with  strong  horizontal 
gradients  of  vertical  velocity  as  we  would  expect. 

Figure  4.1b  shows  the  analyzed  r|  field  and  Figure  4.2b 
shows  this  field  in  the  S-N  crossection.  The  maximum  values 
(positive  and  negative)  for  n  are  found  east  and  west  of  the 
raicroburst  center.  These  maxima  coincide  with  strong  horizon¬ 
tal  gradients  of  vertical  velocity  as  did  the  £  field 


Figure  4.3s  Vertical  velocity  (ra/s)  contoured  on  vector  plot 
of  in-plane  vorticity  field  in  S-N  cross  section  for  1647  MDT, 
14  July  1982. 
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(cross-section  not  shown).  The  actual  centers  are  slightly 
weaker  than  the  £  values  and  are  displaced  in  the  downwind 
(south  through  southeast)  direction  indicating  the  impact  of 
the  environmental  flow. 

The  C  field  is  depicted  in  Figures  4.1c  and  4.2c.  Near 
zero  values  of  C  are  found  below  1  km  at  the  microburst's 
center.  This  would  signify  that  the  flow  is  primarily  diver¬ 
gent  in  this  region,  since  we  have  already  seen  that  £  and  n 

are  also  small  here.  C  increases  above  1  km  to  values  greater 
-3 

than  5  x  10  per  second.  The  microburst's  irrotational  struc¬ 
ture  below  1  km  may  act  to  spin  down  the  positive  C  vorticity 
in  this  area.  This  1  km  level  may  be  a  transition  area  of 
rotation  to  divergence. 

Also  from  Figures  4.1c  and  4.2c,  one  develops  an  overall 
sense  that  this  3-dimensional  vorticity  field  is  largely  a  hor¬ 
izontal  vector  or  at  the  very  least  quasi-horizontal  in  the 
periphery  of  the  raicroburst.  The  values  for  C  are  generally 
one  order  of  magnitude  smaller  than  the  horizontal  components. 
Figure  4.4a  depicts  the  horizontal  vorticity  components  at  0.75 
km  combined  to  form  a  horizontal  vector  field.  A  vortex  ring 
surrounds  the  microburst  indicative  of  the  strong  shear  in  this 
region. 

Fig.  4.4b  illustrates  the  perturbation  pressure  field  from 
Hughes  (1987)  for  the  0.75  km  level.  It  is  plainly  seen  that 
the  vortex  ring  lies  in  a  region  of  low  perturbation  pressure 
which  serves  to  accelerate  the  wind.  Likewise,  the  rotor 


(a) 


w.  Q*  73  KH 
71  " 


16H7  MPT 


P'riiturt  C«W  CS-ND  16H7L  IEj"/* 


Figure  4.4:  Flan  view  of  (a)  vorticity  field  at  0.75  km  and 
(b)  with  the  perturbation  pressure  field  contoured.  S-N  cross 
section  (c)  of  wind  field  with  perturbation  pressure  field  con¬ 
toured  for  1647  MDT,  14  July  1982. 
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itself  is  within  the  low  perturbation  pressure  region  as  viewed 
in  the  S-N  cross  section  in  Fig.  4.4c.  Parsons  efc  al.  (1987), 
Drogemeier  and  Wilhelmson  (1987),  and  Kessinger  et  al.  (1988) 
also  suggest  that  these  vortex  circulations  could  be  partly 
responsible  for  the  observed  high  surface  winds. 

The  vortices  studied  by  Kessinger  et  al.  (1988)  moved 
away  from  the  storm  center.  This  has  lead  them  to  hypothesize 
about  the  possibility  that  variations  in  the  strength  of  the 
downdraft  may  create  separate  centers  of  horizontal  vorticity 
which  then  move  down  and  away  from  the  storm.  The  existence  of 
horizontal  vortices  can  be  inferred  from  Fig.  4.2a  and  Fig. 
4.2b  on  either  side  of  the  microburst.  The  vorticity  centers 
are  more  evident  in  SW-NE  cross  section  (not  shown  for  1647 
MDT,  see  Fig.  4.7  for  this  cross  section  at  1649  MDT) .  Time 
resolution  does  not  allow  tracking  of  these  vorticity  centers, 
but  some  movement  down  and  outward  can  be  inferred  between  1647 
and  1649  MDT  northeast  of  the  microburst.  This  apparent  move¬ 
ment  could  be  a  restablishment  or  a  development  of  vorticity 
lower  and  further  away  as  opposed  to  actual  movement. 

From  the  analyzed  C  field  (Fig.  4.1c),  this  vortex  ring 
oscillates  through  the  plane  at  0.75  km  producing  alternating 
regions  of  positive  and  negative  C  values.  Finally,  the 
strongest  values  of  C  are  located  above  the  1  km  level  in  the 
microburst  area  and  at  the  eastern  edge  of  the  gust  front  while 
still  being  positive  along  the  gust  front,  indicative  of  the 
existence  of  horizontal  shear  in  these  regions. 
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In  order  to  determine  typical  magnitudes  of  the  three  vor- 
ticity  components,  area  means  and  standard  deviations  were  com¬ 
puted  over  the  full  storm  domain  and  the  microburst  domain  (box 
in  Figure  4.1).  These  results  are  listed  in  Table  4.1.  A  com¬ 
parison  of  root  mean  square  values  (standard  deviations)  in  the 
table  shows  that  the  £  and  n  components  are  the  dominant  com¬ 
ponents  at  every  level.  The  £  component  is  slightly  larger 
than  the  n  component  in  most  cases.  The  rotor,  which  effects 
the  vorticity  field  at  levels  0.5  km  through  1  km,  and  the 
predominant  northerly  flow  to  the  northeast  produce  strong 
positive  values  of  £  (Figure  4.1a).  This  results  in  large 
positive  values  for  the  £  areal  means  at  these  levels.  The 
values  for  q  and  C  are  more  varied,  positively  and  negatively, 
and  hence  their  areal  means  are  smaller  than  that  of  £. 

The  same  results  hold  for  the  microburst  domain  with  the 
exception  that  the  n  component  is  strongest  (rms  value)  at  lev¬ 
els  0.75  and  1  km.  Comparing  across  the  domains,  one  finds  the 
maximum  positive  and  maximum  negative  values  of  each  component 
over  the  grid  illustrates  the  significance  of  the  presence  of 

the  microburst,  that  is,  these  values  occur  in  the  microburst 

-2  -2 

domain  predominantly.  Peak  values  of  1.9  x  10  and  -1.6  x  10 

-2  -2 

per  second  for  £  and  1.7  x  10  and  -1.7  x  10  per  second  for 
r|  at  0.75  km  signifies  the  vortex  ring  resides  mainly  at  this 


level  as  was  shown  by  Figure  4.4a. 


Table  4.1:  Area  mean  (<  >)  and  standard  deviation  (rms)  values 
for  each  vorticity  component  over  (a)  storm  domain  and  (b) 
microburst  domain  for  14  July  1982,  1647  MDT.  Units  are 
(  xlO-4)  s  . 

(a) 


Ht(km) 

<i> 

^rras 

<n> 

rms 

<C> 

^rros 

0.25 

-1.8 

45.6 

-3.5 

29.0 

0.87 

15.6 

0.50 

13.5 

46.5 

-1.4 

33.3 

-0.22 

14.9 

0.75 

28.5 

50.5 

-1.8 

47.4 

-1.9 

14.4 

1.00 

19.8 

40.7 

-3.8 

42.5 

-1.9 

17.1 

1.50 

3.5 

41.4 

-2.9 

38.9 

-0.06 

19.9 

2.00 

-6.5 

43.6 

-0.4 

38.1 

1.4 

22.9 

<b) 


Ht(km) 

<S> 

^rras 

<n> 

nrms 

<C> 

^rms 

0  25 

9.9 

46.7 

-3.20 

34.3 

-0.60 

15.9 

0.50 

18.4 

50.5 

-4.60 

47.5 

-3.50 

14.6 

0.75 

24.6 

68.6 

-0.95 

77.1 

-2.70 

13.2 

1.00 

21.2 

55.7 

3.60 

64.0 

-0.45 

19.3 

1.50 

16.5 

53.3 

-0.48 

53.1 

1.20 

24.0 

2.00 

10.7 

53.8 

-7.80 

47.8 

5.00 

31.1 

23 


24 


4.1.1  Budget  Terms;  Simple  Case,  1647  MPT 

Tables  4. 2-4. 7  list  the  results  of  the  computed  values  of 
the  terms  in  each  of  the  component  equations.  They  provide  a 
level  by  level  analysis  of  the  areal  mean  and  standard  devia¬ 
tions  for  each  of  the  four  terms,  horizontal  advection  (HAD), 
vertical  advection  (VAD),  divergence  (DIV),  and  tilting  (TILT) 
for  the  storm  domain  (Tables  4. 2-4.4)  and  the  microburst  domain 
(Tables  4. 5-4. 7).  The  arrows,  <  >,  indicate  area  means. 

Since  the  area  means  for  each  component  of  the  £>  n ,  and  C 
are  at  least  one  order  of  magnitude  smaller  than  their  standard 
deviations  at  each  level,  the  standard  deviations  may  be  con¬ 
sidered  their  typical  values.  Table  4.2  lists  the  budget  terms 
for  while  budget  terms  for  n  and  C  are  found  in  Table  4.3 
and  4.4,  respectively.  HAD  is  the  largest  of  all  the  terms  for 
£  and  q  budgets  with  VAD  and  TILT  being  the  next  higher  terms, 
respectively.  DIV  is  the  smallest  of  all  terms.  TILT  dom¬ 
inates  the  VAD  term  in  the  C  budget  and  has  the  strongest  typi¬ 
cal  value  of  all  terms  within  the  C  budget  at  0.75  km.  It  was 
shown  earlier  that  the  £  component  of  vorticity  was  the  strong¬ 
est  component  and,  therefore,  its  budget  has  the  strongest 
terms  of  all  three  directions. 

The  last  row  in  Tables  4. 2-4. 4  represent  the  vertical 
total  of  the  area  means  for  £  with  HAD,  VAD,  and  TILT  being 
nearly  equal.  HAD  and  TILT  play  opposite  roles  as  source  and 
sink  for  n  (Table  4.3)  as  do  VAD  and  DIV  but  to  a  lesser 
degree.  It  appears  that  the  advection  terms  are  nearly 
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TOTAL  -4.56  3.68  -0.09 
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TOTAL  -0.68  11.60  0.75  -11.20  2.74 


31 


balanced  by  the  divergence  and  tilting  terms  for  this  direc¬ 
tion.  The  final  budget,  C»  indicates  that  VAD  is  the  source 
while  the  other  terms  are  all  acting  to  decrease  positive  C* 
This  follows  from  the  vertical  cross  section  (Figure  4.2c)  in 
that  C  was  positive  above  0.75  km  for  most  of  the  storm  domain. 

Tables  4. 5-4. 7  lists  the  area  mean  and  standard  deviations 
for  the  microburst  domain  budget  terms.  As  expected,  the  typi¬ 
cal  magnitudes  (standard  deviations)  are  larger  in  the  micro¬ 
burst  domain  as  compared  to  those  in  the  storm  domain  for  each 
term  at  virtually  every  level.  The  hierarchy  of  significance 
remains  as  in  the  storm  domain  as  HAD,  VAD,  TILT,  and  DIV  for 
the  £  and  n  budgets.  The  C  budget  again  has  VAD  and  TILT  being 
very  close  in  magnitude. 

The  vertically  totalled  area  means  are  found  in  the  last 
row  of  Table  4. 5-4. 7.  In  Table  4.5,  HAD  acts  to  decrease 
(sink)  S,  while  the  remaining  terms  provide  a  source  of  posi¬ 
tive  £.  VAD  and  TILT  have  the  largest  totals,  but  the  deficit 
attributed  to  HAD  cuts  the  source  due  to  VAD  in  half.  The 
positive  contributions  from  DIV  and  TILT  are  roughly  five  and 
six  times,  respectively,  those  at  the  storm  domain.  This 
points  out  the  significance  of  the  microburst  flow  field 
towards  the  generation  of  £. 

The  n  budget  in  the  microburst  domain  (Table  4.6)  shows 
via  the  vertical  total  row  the  same  trend  as  at  the  storm 
domain.  That  is,  HAD  and  VAD  are  sources  of  while  DIV  and 
TILT  are  sinks.  The  largest  change  occurs  with  VAD  and  DIV 
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between  domains.  VAD  develops  to  become  the  dominant  source 
term  and  DIV  increases  fourfold  to  virtually  match  TILT  as  a 
sink.  Yet  again,  the  sources  remain  due  to  HAD  and  VAD. 

VAD  and  TILT  dominate  the  budget  terms  in  the  C  component 
for  the  microburst  domain.  VAD  is  the  source  and  TILT  is  the 
sink.  The  decrease  in  HAD  seems  to  indicate  the  horizontal 
gradients  of  C  are  weaker  on  average  than  the  vertical  gra¬ 
dients.  Furthermore,  as  VAD  is  acting  to  bring  positive  C 
values  lower  into  the  storm,  TILT  appears  to  single-handedly 
spin  down  the  positive  values.  Figure  4.5  illustrates  the 
opposing  nature  of  these  two  terms  in  profile.  The  result  of 
this  interaction  is  to  keep  the  C  component  of  vorticity  small 
and  hence  the  3-dimensional  vorticity  vector  quasi-horizontal. 

Although  not  evident  from  the  vertical  totals,  HAD  and  VAD 
tend  to  oppose  one  another  in  vertical  profile  plots  of  their 
mean  values.  DIV  and  TILT  then  would  be  required  to  balance 
the  result  of  the  advection  fields  in  order  for  vorticity  to 
remain  unchanged. 

The  last  two  columns  in  Tables  4.2  -  4.7  represent  the 
areal  mean  and  standard  deviations  for  the  respective  tendency, 
i.e.,  d£/dt,  etc.  These  were  computed  from  the  calculated 
values  of  the  four  terms  of  the  vorticity  component  equations 
and  not  measured  directly  so  care  must  taken  in  their  interpre¬ 
tation.  The  tendency  terms  are  less  reliable  than  all  the 
other  terms  in  the  vorticity  budget  equation.  In  addition  to 
the  uncertainities  in  u,  v,  and  w  estimates,  the  solenoidal 


2.0 


(®t)  aaminv 


33 


Fig.  4.5:  Vertical  profile  of  vertical  advection 
(dashed  line)  and  tilting  term  (solid  line)  for  the 
< ;  component  equation  over  tiie  microburst  domain  at 
1647  MDT,  14  July  1982. 
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and  frictional  terms  are  excluded  in  the  budget  calculation. 
So  these  terms  include  the  physical  mechanisms  affecting  the 
respective  vorticity  component  tendency,  but  also,  errors  in 
computing  the  three  wind  components  and  truncation  errors.  The 
truncation  errors  were  minimized  by  using  fourth-order  finite 
difference  schemes.  Lin  et  al.  (1987)  found  errors  in  the 
derived  horizontal  wind  field  to  be  about  1-2  meters  per 
second.  These  errors  are  due  to  the  statistical  uncertainty  of 
the  radial  velocity  estimates  and  geometrical  considerations. 
Estimates  in  vertical  velocity  lead  to  minimal  errors  when  the 
anelastic  continuity  equation  is  integrated  downward  with  a 
variational  adjustment  as  was  down  in  this  storm  (Lin  et  al. , 
1987  and  Lin  and  Coover,  1988).  Further  error  analysis  details 
are  discussed  by  Hughes  (1986),  Lin  et  al.  (1987),  and  Coover 
(1988).  Errors  due  to  the  omission  of  the  solenoidal  and  fric¬ 
tion  terms  in  budget  calculations  will  show  up  in  the  tendency 
terms.  That  is,  the  budget  terms  not  being  balanced  result  in 
a  larger  (positive  or  negative)  tendency  term.  Since  the  ten¬ 
dency  term  was  not  measured  directly,  no  error  can  be  made. 

4.2  Simple  Case;  14  July  1982,  1649  MPT 

Figure  4.6a-4.6c  depict  the  analyzed  vorticity  field  over- 
lain  on  the  wind  field  for  the  0.5  km  level  at  1649  MDT.  Two 
new  gust  fronts  have  developed  (dashed  lines).  £  (Fig.  4.6a) 
has  changed  most  noticeably  along  the  gust  fronts  and  the 
orientation  of  the  field  is  now  NW-SE,  that  is,  parallel  to 
the  fronts,  n  has  also  increased  (Fig.  4.6b)  in  magnitude 
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denoted  (straight  line). 
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mostly  in  the  southwest  region  of  the  microburst  domain.  An 

increase  in  r|  is  also  apparent  by  a  change  in  the  coverage  of 
-3  -1 

the  10  x  10  s  isoline  east  of  the  microburst.  Likewise, 
the  C  component  (Fig.  4.6c)  has  become  more  positive  along  the 
gust  front  to  the  east  of  the  microburst. 

Figure  4.7a-4.7c  is  a  SW-NE  crossection  through  the  center 
of  the  raicroburst  with  £  (Fig.  4.7a),  n  (Fig-  A. 7b),  and  C 
(Fig.  4.7c)  overlain  on  the  wind  field.  The  £  and  n  vorticity 
fields  are  essentially  zero  at  the  microburst  center  now 
through  the  full  depth  and  C  has  become  slightly  positive.  The 
storm  is  in  its  decay  stage  at  this  time  as  evidenced  by  the 
weaker  vertical  velocities  in  the  microburst  region.  The  gust 
front  to  the  east  appears  to  be  better  defined  than  the  western 
gust  front  from  a  look  at  the  winds  in  this  cross  section.  In 
fact,  the  flow  or  horizontal  gradient  of  the  vertical  velocity 
from  the  microburst  center  northeastward  has  produced  a  series 
of  positive  £  vorticity  centers  with  values  in  excess  of  15  x 


The  western  gust  front  is  also  evident  in  this  cross  sec¬ 
tion  by  the  couplet  of  positive  and  negative  f;  centers  1-3  km 
SW  of  the  microburst  center.  Also  of  note  here  is  the  inten¬ 
sity  of  the  horizontal  wind  outflow  in  these  regions.  The  max¬ 
imum  outflow  occurs  below  the  maximum  vorticity  values.  This 
was  also  shown  by  Kessinger  et  al.  (1988)  and  Waranauskas 
(1985)  to  occur  with  their  horizontal  vortex  circulations. 
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As  in  the  1647  MDT  case,  area  means  and  standard  devia¬ 
tions  were  computed  (See  Table  4.8).  £  is  again  the  dominant 
vorticity  component  from  a  standard  deviation  or  root  mean 
square  perspective  in  both  the  storm  and  microburst  domains. 
The  largest  root  mean  square  values  for  £  and  q  now  appear 
lower  in  the  storm,  0.25  km  -  0.75  km,  than  they  were  at  1647 
MDT.  This  is  also  evident  in  Figure  4.7a  and  Figure  4.7b. 
4.2.1  Budget  Terms:  Simple  Case,  1649  MDT 

The  budget  terms  for  the  1649  MDT  time  period  are  listed 
in  terms  of  area  means  and  standard  deviations  over  the  storm 
domain  (Tables  4.9-4.11)  and  for  the  microburst  domain  (Tables 
4.12-4.14).  The  first  item  of  importance  is  the  increased  mag¬ 
nitude  of  the  standard  deviation  for  DIV  in  each  direction  over 
the  1647  MDT  analysis  time  (See  Figures  4.2  and  4.3  for  com¬ 
parison).  These  values  have  increased  at  least  one  order  of 
magnitude  for  each  component  direction  for  levels  below  1  km. 
The  divergent  flow  destroys  positive  C  and,  at  this  time,  the 
vertical  shear  has  become  significant  in  the  lowest  levels  as 
the  storm  continues  to  decay.  The  result  is  that  £  and  q  are 
stronger  in  low  level  while  C  is  strongest  negative  (positive) 
in  low  levels  (higher  levels)  of  the  analysis  volume. 

The  budgets  of  the  three  components  at  1649  MDT  are  simi¬ 
lar  to  those  at  1647  MDT.  The  major  change  is  the  increased 
role  that  DIV  and  TILT  play.  It  seems  the  vorticity  associated 
with  the  1647  MDT  vortex  ring  establishes  itself  closer  to  the 
surface  (See  Table  4.8b).  Figure  4.8  shows  the  £  and  q 


Table  4.8:  Area  mean  (<  >)  and  standard  deviation  (rms)  values 
for  each  vorticity  component  over  (a)  storm  domain  and  (b) 


microburst  domain  for  14  July  1982,  1649  MDT.  Units  are 
(  xlO  )  s 


(a) 


Ht(km) 

<!;> 

rms 

<n> 

n  „ 
rms 

<c> 

^rms 

0 . 25 

15.10 

78.5 

-12.90 

47.8 

-2.80 

28.7 

0.50 

26.80 

79.7 

-16.40 

54.8 

-1.30 

20.9 

0.75 

35.50 

70.0 

-14.50 

52.6 

-0.26 

20.3 

1.00 

15.50 

47.3 

-17.90 

36.1 

-0.22 

23.4 

1.50 

0.62 

42.1 

0.94 

28.7 

1.10 

23.7 

2.00 

-5.60 

47.6 

0.18 

29.1 

2.40 

22.5 

Ht(km) 

<*;> 

^rms 

<n> 

n 

rms 

<C> 

^rms 

0.25 

-10.20 

92.1 

12.2 

64.4 

6.8 

30.8 

0.50 

-0.41 

88.5 

12.6 

70.0 

2.4 

20.2 

0.75 

23.40 

71.8 

7.6 

63.5 

-1.2 

19.2 

1.00 

27.50 

50.0 

1.4 

47.8 

3.6 

19.6 

1.50 

14.20 

37.2 

-0.5 

33.7 

6.1 

24.1 

2.00 

15.70 

41.0 

-1.6 

26.6 

6.9 

2.9 
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divergence  profiles  for  the  microburst  domain  at  both  analysis 
times.  The  tilting  profiles  show  a  similar  trend  and,  hence, 
are  not  shown  here.  The  TILT  and  VAD  terms  are  again  compen¬ 
sating  for  the  C  budget  so  the  vorticity  remains  quasi¬ 
horizontal  . 


4.3  Complex  Case;  5  August  1982 

The  analyzed  vorticity  field  for  1845  MDT  at  the  0.5  km 
level  is  presented  on  the  wind  field  in  Figure  4.9.  This 
storm  is  more  complicated  than  the  previous  case  of  14  July 
1982.  As  depicted  in  Fig.  4.9,  this  storm  produced  a  micro¬ 
burst,  Ml  (-2, -24), a  gust  front,  northeast  through  southeast 
(dashed  line)  of  Ml,  and  an  enhanced  downdraft,  D  (-6,-19). 
The  complex  interactions  of  the  flow  fields  of  the  enhanced 
downdraft  and  Ml  combine  to  create  a  large  elongated  region  of 
the  two  horizontal  components  of  vorticity  from  (-7,-17)  to  (- 
2,-21).  This  area  corresponds  to  an  area  of  low  perturbation 
pressure  gradient  as  retrieved  by  Coover  (1988)  (not  shown), 
and  is  a  feature  found  in  horizontal  vortex  circulations  (Kess- 
inger  et  al.,  1988).  A  circulation  is  not  found  here,  but 
rather,  the  horizontal  vorticity  is  high.  The  actual  shear 
then  is  not  apparent  simply  by  examining  the  Doppler  derived 
wind  field.  The  raicroburst  domain  (small  box)  encompasses  Ml, 
a  large  portion  of  the  gust  front,  and  the  misocyclone,  c^. 
The  inflowing  environmental  air  is  evident  east  of  Ml  at  the 
bulge  of  the  gust  front. 

4.3.1  X-direction  Vorticity  analysis  (£) 

£  (Fig.  4.9)  is  overwhelmingly  negative  in  the  microburst 
domain  at  all  levels  with  a  pattern  as  shown  here  for  the  0.5 
km  level.  £  is  positive  to  the  northwest  as  depicted  in  the 
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NW-SE  crossection  (A-B,  Fig.  4.9)  in  Figure  4.10.  Above  0.5 
km,  £  is  positive  in  this  northwest  quadrant.  It  is  in  this 
region  that  the  divergent  outflow  of  Ml  is  strongest  in  the 
lowest  levels. 

The  reason  for  the  distribution  of  £  as  largely  negative 
south  of  M^  and  positive  north  can  be  seen  from  Fig.  4.10.  The 
inflowing  air  velocity  increases  with  height  south  of  M^ ,  i.e., 
3v/3z  is  positive  and  decreases  strongly  in  the  north.  This 
results  in  the  given  1845  MDT  £  distribution.  By  1850  MDT,  the 
only  real  change  is  the  most  shallow  layer  of  maximum  outflow 
winds  to  the  north,  i.e.,  3v/3z  is  large  negative  and  there¬ 
fore,  £  is  large  positive  (3w/3y  is  weak). 

£  values  change  very  little  between  1845  and  1847  MDT 
(Fig.  4.11).  From  Tables  4.15  and  4.16  only  a  slight  increase 
in  the  mean  values  of  £  in  low  levels  below  1  km  over  both 
domains  and  a  decrease  above  1  km  is  detected.  Root  mean 
square  values  remain  the  same. 

By  1850  MDT,  £  continues  virtually  unchanged  at  0.5  km 
(Fig.  4.12)  except  for  the  northwest  quadrant  of  the  microburst 
domain.  Over  the  storm  domain,  the  elongated  region  of  strong 
£  values  mentioned  at  1845  MDT  from  north  of  D,  to  north  of  Ml 
has  intensified.  The  perturbation  pressure  field  has  also 
lowered  in  the  same  regions  (Fig.  4.13).  A  band  of  strong  posi¬ 
tive  £  values  10  x  10  ^  s  *  to  20  x  10  ^  s  *  to  15  x  10  ^  s  * 
lies  west  southwest  to  north  of  Ml.  The  wind  field  in  this 
region  curves  anticyclonically  (at  0.5  km)  emanating  from  the 


ALTITUDE  CKM)  ALTITUDE  CKM3 


(b)  XI  CNV-SE3  1850L 


10  m/  s 


Fig.  4.10:  HW-SE  vertical  cross  section  of  £  vorticity  con¬ 
toured  on  wind  field  for  (a)  1845  MDT  and  (b)  1850  MDT,  5 
August  1982.  Location  of  microburst,  Ml  is  indicated. 


KM  EAST  OF  CP-2 


Fig.  4.11:  Same  as  Fig.  4.9  except  for  1847  MDT,  5 
August  1982. 


52 


Table  4.16:  Same  as  Table  4.15  except  for  the  microburst  domain 


(a) 

Ht(km) 

<£> 

^rms 

<n> 

^  rms 

<C> 

^rms 

0.25 

-6.50 

6.5 

-2.50 

11.6 

0.46 

4.7 

0.50 

-4.30 

6.4 

-1.10 

10.0 

0.98 

4.4 

0.75 

-1.00 

7.3 

0.47 

5.8 

1.2 

3.9 

1.00 

0.78 

6.5 

0.74 

5.1 

1.2 

4.2 

1.25 

1.40 

8.7 

0.88 

5.8 

1.0 

4.2 

(b) 

Ht(km) 

<o 

rms 

<n> 

11  rms 

<C> 

^rms 

0.25 

-4.50 

6.1 

-3.3 

9.4 

0.49 

4.9 

0.50 

-3.80 

6.1 

-3.2 

9.4 

1.17 

4.8 

0.75 

-1.60 

7.7 

-1.7 

8.4 

1.34 

4.6 

1.00 

0.08 

7.7 

1.3 

4.6 

1.26 

4.4 

1.25 

0.43 

8.9 

-1.0 

2.9 

1.59 

4.7 

(c) 

Ht(km) 

<£> 

^rras 

<n> 

^rms 

<c> 

^rms 

0.25 

-1.87 

9.2 

-3.62 

9.6 

0.48 

4.1 

0.50 

-0.82 

8.1 

-2.47 

8.7 

0.99 

4.4 

0.75 

0.70 

8.2 

-0.02 

6.9 

1.22 

4.6 

1.00 

2.22 

7.3 

2.15 

5.0 

1.40 

4.0 

1.25 

3.19 

7.6 

3.06 

5.0 

1.48 

3.8 
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high  perturbation  pressure  dome  (0.2  mb)  associated  with  Ml 
into  a  lower  pressure  region  between  both  Ml  and  the  high  pres¬ 
sure  dome  associated  with  D. 

Table  4.15c  and  4.16c  also  shows  that  at  1850  MDT,  £ 
increases  in  the  area  mean  at  all  levels  in  the  storm  domain 
and  microburst  domain.  In  fact,  the  largely  negative  5  mean 
values  in  low  levels  that  existed  at  1845  and  1847  MDT  have  all 
but  been  reversed  in  the  mean  by  1850  MDT.  This  indicates  the 
strengthening  of  a  positive  £  field  as  shown  in  Figure  4.12 
northwest  of  Ml.  Figure  4.10  illustrates  the  stronger  vertical 
shear  at  1850  MDT  responsible  for  the  stronger  positive  £ 
values . 

4.3.2  £  Budget  Terms 

Figure  4.14  illustrates  the  change  in  <£>  (where  <  > 
represents  area  mean)  through  each  time  period  in  vertical  pro¬ 
file  for  the  microburst  domain.  It  can  be  seen  that  £  is 
increased  each  time  below  0.5  km  and  consistently  increases 
with  height  at  each  time  period.  Above  0.5  km,  <£>  at  1845  MDT 
decreases  slightly  by  1847  MDT,  but  is  again  increased  at  1850 
MDT. 

Table  4.17  lists  the  area  mean  and  standard  deviation 
values  for  each  of  the  four  budget  terms.  At  1845  MDT  (Table 
4.17a)  in  the  storm  domain,  a  balance  must  be  accomplished 
between  HAD  and  the  combined  effect  of  DIV  and  TILT  since  VAD 
is  small.  A  balance  is  not  met  and  by  comparing  vertical 


0.75 


(mi)  aaniixiv 


al  profile  of  area  mean  values  of 
1845  (dashed  line),  1847  (dotted 
iT  (solid  line)  5  August  1982. 


Table  A. 17:  Area  mean  (<  >)  and  standard  deviation  (o)  values 
for  the  terms  of  the  x-component  vorticity  equation  over  the 
storm  domain  for  (a)  1845,  £b)  18A7,  and  (c)  1850  MDT,  5  Aug 

1982.  Units  are  (  xlO  )  s 


(a) 

Ht(km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

a 

<DIV> 

DIV 

0 

<TILT> 

TILT 

0 

0.25 

4.20 

39.5 

0.009 

2.7 

-1.30 

16.9 

1.20 

19.0 

0.50 

3.30 

34.8 

0.33 

8.7 

-1.30 

9.9 

1.70 

17.3 

0.75 

0.92 

30.6 

0.59 

10.0 

-0.63 

9.5 

-0.52 

16.8 

1.00 

-0.61 

30.1 

0.03 

13.0 

0.49 

10.1 

-1.20 

16.2 

1.25 

-3.50 

44.3 

-0.42 

15.9 

0.72 

13.7 

-1.70 

20.9 

TOTAL 

4.31 

0.54 

-2.02 

-0.52 

(b) 

Ht(km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

0 

<DIV> 

DIV 

o 

<TILT> 

TILT 

o 

0.25 

4.3 

45.1 

0.15 

2.6 

-1.70 

14.6 

2.3 

20.3 

0.50 

2.8 

36.8 

0.53 

8.9 

-1.30 

12.1 

2.2 

17.5 

0.75 

0.6 

33.8 

0.63 

11.0 

-0.94 

11.4 

1.2 

17.1 

1.00 

1.2 

29.9 

0.46 

15.5 

0.21 

10.2 

0.15 

10.5 

1.25 

1.6 

35.5 

0.51 

19.4 

-0.94 

11.3 

0.94 

18.9 

TOTAL 

10.5 

2.28 

-4.67 

6.79 

(c) 

Ht(kra) 

<HAD> 

HAD 

a 

<VAD> 

VAD 

0 

<DIV> 

DIV 

0 

<TILT> 

TILT 

o 

0.25 

5.0 

55.0 

0.04 

2.8 

-2.40 

20.8 

4.65 

24.4 

0.50 

2.3 

39.2 

-0.11 

9.5 

-1.20 

12.0 

2.98 

21.7 

0.75 

-1.2 

31.8 

-0.06 

12.2 

0.27 

9.4 

-0.65 

20.5 

1.00 

-1.6 

26.1 

1.04 

12.6 

1.40 

12.0 

-0.95 

18.3 

1.25 

-2.7 

29.9 

1.62 

13.4 

1.63 

10.2 

-1.60 

20.3 

TOTAL 

1.8 

2.53 

-0.30 

3.83 
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totals,  a  net  Increase  of  £  is  suggested.  By  1847  MOT  (Table 
4.17b),  the  contribution  from  VAO  has  increased  slightly  (level 
by  level)  but  is  still  one  order  smaller  in  the  mean  than  the 
remaining  terms.  At  this  time,  HAD  and  TILT  acting  positively 
(source)  at  low  levels  overwhelms  the  loss  through  DIV.  Across 
each  time  period,  it  is  apparent  that  the  DIV  term  is  the  sole 
generator  of  negative  £.  The  remaining  terms  generate  positive 
S  with  the  strongest  generation  occuring  below  0.75  km. 

In  the  smaller  microburst  domain  (Table  4.18),  the  area 
mean  and  root  mean  square  (standard  deviation)  values  are 
larger  than  at  the  storm  domain.  From  Fig  4.9  it  is  evident 
that  due  to  the  overwhelming  negative  £  field,  HAD  would  act  as 
a  sink  as  the  table  indicates.  VAD  counters  by  acting  as  a 
source.  TILT  is  a  source  at  each  time  period  with  maximum  con¬ 
tribution  ocurring  at  1847  MDT.  The  large  horizontal  gradients 
of  £  create  a  very  strong  HAD  terra  leading  to  removal  of  posi¬ 
tive  E  within  the  microburst  domain  over  time.  VAD  and  TILT 
counteract  this  removal  mechanism  and  dominate  at  1845  MDT  and 
1847  MDT.  By  1850  MDT  when  the  £  gradient  is  strongest,  DIV 
acts  with  VAD  and  TILT  to  attempt  a  balance.  As  the  £  field 
gets  stronger,  the  DIV  term  especially  to  the  north  of  M^ 
becomes  a  strong  source. 

4.3.3  Y-direction  Vortlcity  Analysis  (n) 

Between  1845  and  1850  MDT,  <n>  values  remain  predominantly 
negative  (Table  4.15  and  4.16)  over  the  storm  domain  and  below 
1  km  in  the  raicroburst  domain.  A  tight  gradient  of  r| 


Table  4.18:  Same  as  Table  4.17  except  for  microburst  domain 


(a) 


Ht(kis) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

0 

<DIV> 

D1V 

0 

<TILT> 

TILT 

0 

0.25 

-9.40 

45.8 

1.00 

5.3 

4.5 

33.4 

-5.90 

34.6 

0.50 

-4.80 

36.9 

3.70 

16.9 

-3.8 

15.5 

7.30 

28.7 

0.75 

-4.70 

46.1 

1.90 

17.3 

-1.4 

11.8 

0.16 

18.9 

1.00 

0.46 

37.5 

-0.16 

24.2 

-2.0 

17.1 

3.20 

25.5 

1.25 

2.00 

53.4 

0.68 

31.8 

-1.7 

19.5 

6.20 

29.6 

TOTAL 

-16.40 

7.12 

-4.4 

10.96 

(b) 

Ht(km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

0 

<DIV> 

DIV 

0 

<TILT> 

TILT 

a 

0.25 

-6.40 

66.7 

0.46 

5.1 

-0.39 

26.1 

-1.70 

32.7 

0.50 

-11.00 

58.8 

2.49 

15.9 

-1.79 

18.5 

0.43 

26.5 

0.75 

-11.60 

56.9 

6.02 

16.3 

-2.16 

15.7 

3.22 

22.4 

1.00 

-4.00 

56.1 

4.94 

20.0 

-1.14 

17.6 

5.61 

27.7 

1.25 

0.96 

67.4 

2.35 

25.2 

-3.39 

24.4 

8.85 

38.2 

TOTAL 

-32.04 

16.26 

-8.87 

16.41 

(c) 

Ht(km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

a 

<DIV> 

DIV 

a 

<TILT> 

TILT 

a 

0.25 

-16.1 

100.9 

2.4 

5.6 

2.8 

42.0 

0.25 

41.2 

0.50 

-20.9 

66.5 

7.7 

18.7 

1.3 

18.7 

4.55 

34.5 

0.75 

-19.0 

49.7 

7.7 

22.2 

2.3 

13.6 

3.51 

31.9 

1.00 

-10.1 

36.6 

4.8 

20.7 

3.6 

18.1 

3.11 

22.3 

1.25 

-6.1 

37.9 

4.7 

19.0 

3.3 

13.9 

3.31 

18.1 

TOTAL 

-72.2 

27.3 

13.3 

14.73 
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results  from  the  strong  outflow  of  Ml  eastward  at  low  levels 
and  the  strong  inflowing  air  especially  east  of  Ml.  This  helps 
locate  the  gust  front  associated  with  Ml  at  1845  and  1847  MDT 
as  well  as  a  new  gust  front  that  develops  at  1850  MDT  associ¬ 
ated  with  the  M3  (Figures  4.15,  4.16,  and  4.17). 

The  most  interesting  changes  for  n  occur  within  the  micro¬ 
burst  domain  as  maximum  positive  and  negative  n  centers  pivot 
around  Ml  and  become  oriented  northwest  to  southeast  through 
the  center  of  Ml  with  positive  n  to  the  northwest  (Fig.  4.17). 
The  tight  horizontal  gradient  of  n  values  along  the  gust  front 
weakens  slightly  with  time  and  it  appears  in  crossection  (Fig. 
4.18)  that  the  slope  of  the  front  is  most  shallow  at  1850  MDT 
in  fact  approaching  horizontal  above  and  southeast  of  Ml.  This 
is  to  say  that  negative  n  is  being  generated  and  is  being 
spread  out  to  the  southeast  in  the  lower  levels  of  the  storm. 
Tables  4.15  and  4.16  shows,  especially  in  the  microburst 
domain,  that  positive  n  is  being  generated  above  0.75  km  and 
destroyed  below  0.75  km. 

The  n  field  in  Figure  4.15  has  its  strongest  negative 
-3  -1 

value,  -32  x  10  s  ,  at  this  level  (0.5  km)  just  east  of  Ml. 
The  gust  front  is  closest  to  Ml  at  this  point.  The  n  field 
parallels  if  not  defines  the  location  of  the  gust  front  with 
its  very  tight  gradient  adjacent  to  the  gust  front. 


Fig.  4.15:  Same  as  Fig.  4.9  except  for  n  vorticity 
field. 
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Fig.  4.16:  Same  as  Fig.  4.11  except  for  n  vorti- 
city. 
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Fig.  4.17:  Same  as  Fig.  4.12  except  for  n  vorti- 
city. 
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ALTITUDE  CKMI  ALTITUDE  CKM3 


(a)  ETA  CNV-SE3  18H5L 


10  m/s 


Figure  4.18:  NW-SE  vertical  cross  section  of  n  vorticity  on 
wind  field  for  (a)  1845  and  (b)  1850  MDT,  5  August  1982.  Loca¬ 
tion  of  Ml  is  indicated. 
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4.3.4  n  Budget 

Figure  4.19  displays  the  vertical  profile  of  <n>  at  the 
microburst  domain.  This  change  is  more  significant  than  at  the 
storm  domain.  <n>  experiences  a  decrease  at  the  1847  MDT  time 
period.  <n>  increases  again  by  1850  MDT  but  does  not  exceed 
its  initial  1847  MDT  vertically  averaged  value.  The  tilting 
term  seems  to  play  the  biggest  role  in  making  this  decrease 
occur.  Table  4.19  shows  the  budget  terms  for  n-  HAD  and  VAD 
and  likewise,  DIV  and  TILT  act,  in  general,  to  oppose  each 
other. 

4.3.5  2-direction  Vorticity  Analysis  (Q 

Figure  4.20  shows  the  C  field  in  plan  view  at  0.5  km,  and 
Tables  4.15  and  4.16  list  the  area  means  standard  deviations  by 
level  for  both  domains.  From  Table  4.15,  it  is  apparent  that  C 
changes  very  little  statistically  at  the  storm  domain.  The 
largest  changes  occur  in  the  microburst  domain  (Fig.  4.7). 

From  Fig.  4.20,  two  areas  of  positive  C  are  evident.  One 

-3  -1 

area  west  of  Ml  (area  >  4  x  10  s  )  at  0.5  km  and  the  other 

-3  -1 

along  the  gust  front  with  a  maximum  (14  x  10  s  )  associated 
with  the  mesocyclonic-1 ike  vortex,  c^,  located  at  the  southern 
end  of  this  gust  front  at  1845  MDT. 

By  1847  MDT  (Fig.  4.21),  the  area  along  the  gust  front  has 
been  split  as  C  is  diminished  just  south  of  Ml  and  this  contin¬ 
ues  at  1850  MDT  (Fig.  4.22).  The  raesocyclone-1 ike  vortex,  c^, 
weakens  with  time. 


Table  4.19:  Area  mean  (<  >)  and  standard  deviation  (o)  values 
for  the  terms  of  the  y-component  vorticity  equation  over  the 
microburst  domain  for  ,(a)  1845,  (b)  1847,  (c)  1850  MDT,  5  Aug 
1982.  Units  are  (  xlO  )  s 


(a) 


Ht(km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

a 

<DIV> 

DIV 

a 

<TILT> 

TILT 

a 

0.25 

11.50 

56.1 

-0.19 

4.2 

-0.32 

48.2 

-0.23 

40.1 

0.50 

13.90 

60.5 

-1.90 

16.9 

-2.20 

14.1 

1.90 

30.5 

0.75 

9.00 

53.9 

-4.30 

28.3 

-1.80 

10.2 

3.60 

24.4 

1.00 

0.82 

37.6 

-2.10 

23.4 

-0.80 

12.1 

3.90 

22.3 

1.25 

-4.20 

37.4 

0.64 

16.6 

0.22 

9.8 

0.90 

24.3 

TOTAL 

31.02 

-7.85 

-4.90 

10.07 

(b) 

Ht(kra) 

<HAD> 

HAD 

a 

<VAD> 

VAD 

0 

<DIV> 

DIV 

a 

<TILT> 

TILT 

a 

0.25 

0.79 

50.1 

-0.25 

5.7 

3.5 

39.2 

-16.5 

40.8 

0.50 

6.64 

59.5 

-0.02 

18.5 

1.3 

27.4 

-6.5 

29.9 

0.75 

5.23 

62.7 

3.36 

24.0 

-1.8 

21.0 

3.7 

26.9 

1.00 

3.75 

51.6 

2.71 

26.9 

-1.7 

10.4 

8.8 

23.8 

1.25 

1.95 

35.7 

-0.33 

27.1 

-1.4 

8.4 

8.0 

32.0 

TOTAL 

18.36 

5.47 

-0.1 

-2.5 

Cc) 

Ht(kra) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

0 

<D1V> 

DIV 

0 

<TILT> 

TILT 

a 

0.25 

-14.30 

46.5 

2.00 

5.2 

5.10 

35.7 

-8.0 

38.2 

0.50 

-9.00 

44.6 

6.20 

18.5 

5.70 

22.5 

-7.0 

31.2 

0.75 

-0.41 

29.8 

4.70 

29.8 

-0.33 

11.9 

4.1 

28.5 

1.00 

-0.34 

26.9 

-0.61 

26.9 

-0.25 

10.3 

5.3 

22.3 

1.25 

-3.10 

19.2 

-0.71 

19.2 

1.89 

9.3 

3.1 

20.7 

TOTAL 

-27.15 

11.58 

12.11 

-2.5 
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Figure  4.20:  Same  as  Fig.  4.9  except  C  field  and  W-E  cross 
section  is  depicted. 
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4.21:  Same  as  Fig.  4.20  except  for  1847  MDT,  5  August 
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-3  -1 

The  area  of  C  greater  than  12  x  10  s  tracks  northwestward 
with  time,  weakens  as  it  approaches  Ml,  and  merges  with  the 
positive  C  area  to  the  west  of  Ml. 

A  west-east  cross  section  (Fig.  4.23)  shows  that  initially 
the  area  west  of  Ml  is  shallow  extending  only  up  to  0.5  km.  By 
1847  MDT  and  continuing  at  1850  MDT,  this  area  experiences  a 
spin  up  of  vorticity  from  the  top  down.  A  new  misocyclone,  c2, 
develops  here  by  1850  MDT  with  maximum  C  vorticity  occuring  at 
1.25  km  of  14  x  10~3  s'1. 

-3  -1 

C  is  a  maximum,  14  x  10  s  ,  in  the  mesocyclone-like 
vortex  and  is  positive  along  the  gust  front.  Another  region  of 
postive  C  values  is  located  west  of  Ml.  The  W-E  cross  section 
through  the  microburst  domain  in  Figure  4.23a  indicates  the 
positive  area  west  of  Ml  is  a  low  level  feature  while  the  area 
along  the  gust  front  extends  through  several  levels  and  slopes 
westward  with  height,  i.e.,  with  the  updraft  in  this  area.  In 
the  column  directly  above  Ml,  C  is  negative  below  0.75  km  and 
positive  above  although  the  values  are  relatively  small,  i.e., 
near  zero  throughout.  This  signifies  that  the  flow  is  largely 
irrotational  through  the  horizontal  plane. 

Table  4.15  lists  the  area  mean  and  standard  deviation 
values  for  each  component  of  vorticity  by  level  for  the  storm 
domain.  By  comparing  area  mean  values  of  £,  n»  and  C>  it  is 
seen  that  £  and  £  are  the  predominantly  positive  over  the  storm 
domain  with  the  strongest  values  above  0.5  km.  p  values  have 
negative  values  at  the  lower  levels,  0.25  km  and  0.5  km. 


(a)  ZETA  CV-ED  18^5L 


10  m/s 


(b)  ZETA  CV-ED  18^7L 


Y  =  -25  -25  -25 

Figure  4.23:  W-E  cross  section  of  C  vorticity  at  (a)  1845,  (b) 
1847,  and  (c)  1850  MDT,  5  August  1982.  Gust  front  is  indicated 
by  dashed  line. 
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In  the  raicroburst  domain  (Table  4.16),  £  and  n  are  predom¬ 
inantly  negative  in  the  mean  below  0.75  km  (O  and  0.5  km  (n) 
while  C  continues  to  be  postive  throughout  with  the  larger  mean 
values  at  the  higher  levels  of  the  microburst  domain. 

4.3.6  C  Budget 

The  C  budget  terms  in  the  storm  domain  for  virtually  all 
times  were  found  to  be  less  than  those  in  the  microburst 

domain.  Area  mean  values  were  typically  on  the  order  of  1  x 

7  2  _  ^  _  2 
10  s  in  the  storm  domain  and  1  x  10  s  in  the  microburst 

domain  indicating  the  many  more  complex  and  opposing  interac¬ 
tions  at  the  storm  domain.  Typical  values,  while  being  of  the 
same  order  in  both  domains  were  larger  in  the  microburst  domain 
and  extended  through  more  depth. 

The  microburst  domain  £  budget  terms  are  presented  in 
Table  4.20  for  all  analysis  times.  The  last  row  is  a  vertical 
sum  of  the  area  means  of  each  term.  This  provides  a  look  at 
which  terms  act  as  a  source  or  sink  for  £.  It  can  be  seen  that 
DIV  and  TILT  become  important,  especially  at  1847  and  1850  MDT 
at  0.75  to  1.25  km.  Both  of  these  terms  act  positively  to  gen¬ 
erate  C-  So  it  appears  that  the  contribution  from  DIV  is  crit¬ 
ical  to  the  formation  and  development  of  the  misocyclone. 

Figure  4.24  shows  a  vertical  profile  of  area  means  for  C 
at  each  analysis  time  for  the  microburst  domain.  The  <C>  field 
experiences  virtually  no  change  with  time  while  maintaining  a 
net  positive  component  at  all  levels. 


Table  4.20 

Are  a 

mean 

(<  >) 

and  standard 

deviation  (a) 

values 

of  the  terms  of 

the  z 

-component 

vorticity  equation  over  the 

microburst 
Aug  1982. 

domain  for 
Units  are  ( 

(a). 1845 
xlO"6)  s'2. 

(b)  1847 

,  and 

(c)  1850 

MDT,  5 

(a) 

Ht (km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

0 

<DIV> 

DIV 

a 

<TILT> 

TILT 

O 

0.25 

-1.80 

28.6 

0.23 

3.4 

4.90 

14.4 

1.80 

9.4 

0.50 

2.80 

30.3 

1.5 

8.5 

-0.03 

7.9 

-0.85 

11.6 

0.75 

1.90 

26.9 

1.99 

15.2 

-0.30 

8.4 

-2.10 

13.9 

1.00 

0.54 

28.6 

2.1 

11.2 

0.36 

10.1 

-1.30 

15.4 

1.25 

-1.90 

24.1 

2.1 

22.7 

1.70 

7.0 

-0.83 

24.8 

TOTAL 

1.54 

7.92 

6.63 

-3.28 

(b) 

Ht(km) 

<HAD> 

HAD 

0 

<VAD> 

VAD 

0 

<DIV> 

DIV 

0 

<TILT> 

TILT 

0 

0.25 

4.90 

36.8 

-0.  13 

3.9 

-0.75 

16.5 

-1.9 

12.6 

0.50 

5.70 

27.0 

-0.58 

9.4 

-2.26 

9.1 

-1.6 

16.4 

0.75 

0.92 

30.9 

-2.61 

13.9 

0.16 

9.5 

1.9 

16.6 

1.00 

-4.75 

27.9 

-1.58 

10.5 

2.67 

10.3 

4.5 

13.3 

1.25 

-6.50 

26.1 

0.45 

14.4 

4.00 

9.3 

4.6 

16.2 

TOTAL 

0.27 

-4,45 

3.82 

7.5 

(c) 

Ht(km) 

CHAO 

HAD 

0 

<VAD> 

VAD 

0 

<DIV> 

DIV 

o 

<TILT> 

TILT 

a 

0.25 

8.39 

34.5 

-0.11 

4.89 

-4.40 
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Fig.  A. 24:  Same  as  Fig.  4.14  except  C  vorticlty  is 
depicted . 


5.  Conclusions 


Horizontal  vorticity  maxima/miniraa  were  found  in  regions 
of  strong  horizontal  gradients  of  vertical  velocity  above  the 
lowest  layers.  In  the  lower  layers,  but  outside  the  immediate 
center  of  the  microburst,  these  horizontal  vorticity  centers 
were  enhanced  by  the  vertical  shear  of  the  horizontal  winds. 

Vertical  vorticity,  C>  was  positive  above  about  1  km  for 
the  14  July  1982,  simple  case  and  virtually  zero  or  negative 
below.  The  microburst  flow,  then,  appeared  to  be  irrotational . 

Similar  results  for  the  complex  case  were  found.  The  5 
August  case  had  a  stronger  environmental  flow.  Its  microburst, 
,  was  less  downward  directed,  but  rather,  directed  outward 
(to  the  north)  and  down.  The  horizontal  gradient  of  vertical 
velocity  suffered  in  strength  due  to  these  features.  The  hor¬ 
izontal  components  are  maximized  along  the  gust  front  where 
convergence  is  important  and  when  begins  to  decay  and  its 
outflow  is  teamed  with  that  of  . 

The  simple  case  of  14  July  1982  illustrated  the  existence 
of  a  vortex  ring,  envisioned  by  Fujita  (1985),  surrounding  the 
microburst  downflow.  This  ring  descended  from  0.75  km  at  1647 
MDT  to  0.5  km  at  1649  MDT.  The  effect  was  to  create  increased 
surface  winds  SW  and  NE  (along  the  gust  fronts)  as  the  ring 
descended.  The  increased  surface  winds  and  the  development  or 
reformation  of  the  gust  front  at  1649  MDT  led  to  the  increased 


78 


79 


significance  of  the  divergence  term  in  creating  positive  and 
negative  valued  horizontal  vorticity  areas.  That  is,  the 
increased  surface  winds  create  a  larger  vertical  wind  shear 
which  implies  larger  horizontal  vorticity  components. 

The  vortex  ring  and  rotor  does  have  low  perturbation  pres¬ 
sure  associated  with  it  as  Kessinger  et  al.  (1988)  found  with 
their  actual  pressure  measurements  via  a  PAM  and  their  pressure 
perturbation  analysis.  The  low  pressure  has  been  deemed 
responsible  for  accelerating  the  surface  winds  as  has  been  men¬ 
tioned.  The  derived  perturbation  pressure  field  describes  the 
necessary  pressure  gradient  to  account  for  the  wind  field.  The 
existence  of  the  vortex  ring  explains  the  location  of  the  low 
pressure  and  maximum  surface  winds.  Finally,  it  was  shown  that 
the  overall  vorticity  vector  was  quasi-horizontal  in  this  case. 
The  z-direction  vorticity  budget  illustrated  the  opposing 
nature  of  the  vertical  advection  and  tilting  terms.  Hence,  C 
was  destroyed  through  the  tilting  term  as  vertical  advection 
attempted  to  build  C  from  the  top  down. 

The  complex  case  of  5  August  1982  showed  the  evolution  of 
the  misocyclone.  It  was  seen  that  the  positive  C  was  generated 
from  the  top  down  as  Kessinger  et  al.  (1988)  and  Ray  et  al. 
(1975)  also  suggr  ed.  The  lack  of  sufficient  convergence  in 
the  low-levels  of  the  storm  perhaps  destroyed  the  opportunity 
for  the  first  circulation  along  the  gust  front  to  spin  up. 
Another  possible  explanation  is  that  M^'s  diverging  flow  served 
to  lessen  the  existing  positive  vertical  vorticity,  C»  in  low- 
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levels  as  witnessed  by  the  disappearance  of  the  raesocyclone- 
like  vortex,  c^.  As  the  diverging  flow  continued,  c^  weakened 
and  disappeared.  The  opposing  flow  of  M^  towards  the  south  is 
not  strong  and  does  not  allow  sufficient  convergence. 

Further,  C  (positive  values)  sloped  westward  over  the 
microburst  with  height  and  allowed  the  VAD  terra  to  act  as  a 
source.  While  the  TILT  term  acted  as  a  sink  at  1845  MOT,  it 
began  to  generate  C  at  the  higher  levels  at  1847  and  1850  MDT. 
The  final  contribution  from  the  DIV  term,  acting  strongly  posi¬ 
tive  particularly  at  1847  and  1850  MDT,  allowed  the  spin  up  of 
C  to  continue  with  the  misocyclone,  resulting. 

Vorticity  budget  analyses  were  also  conducted  for  each 
direction  for  each  case.  The  magnitude  of  the  budget  terms 
over  the  raicroburst  domain  have  been  shown  to  be  larger  than 
those  over  the  larger  storm  domain.  Hence,  the  existence  of  the 
microburst  enhances  the  magnitude  of  the  vorticity  components 
in  the  near  raicroburst  region  of  the  storm.  This  was  true  for 
both  cases  and  for  each  component  of  vorticity. 

Lastly,  some  measure  of  comparison  between  the  cases 
should  be  mentioned.  All  three  components  of  vorticity  are 
stronger  in  a  root  mean  square  comparison  in  the  5  August  case 
at  the  raicroburst  domain.  The  stronger  inflowing  environmental 
air  allowed  for  greater  vertical  shear  and  larger  horizontal 
vorticity  components. 

In  a  budget  sense  and  since  the  components  were  larger  in 
the  5  August  case,  the  budget  terms  should  be  larger  for  the  5 
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August  case  than  for  the  14  July  case.  This  was  found  to  be 
true  with  special  notice  to  the  considerably  larger  DIV  term  in 
the  5  August  case,  especially  in  the  z-direction.  It  is 
believed  that  this  term  contributed  greatly  to  the  development 
of  the  misocyclone. 

There  were  no  apparent  consistent  sources  or  sinks  between 
cases.  That  is,  for  example,  if  the  tilting  term  was  a  source 
for  one  component  at  one  time,  it  was  not  necessarily  a  source 
at  a  later  time  or  for  a  subsequent  case. 

Further  study  is  needed  to  determine  why  the  rotors 
develop  in  some  cases  and  not  in  others.  It  may  be  a  combina¬ 
tion  of  convergence  and  opposing  flow.  That  is,  for  the  14 

July  case,  the  rotor  did  not  form  south  of  the  microburst  where 

-3  -1 

the  divergence  magnitude  (-12  x  10  s  )  is  the  same  as  in  the 
north  where  the  rotor  did  form.  Therefore,  convergence  alone 
is  not  enough,  but  the  flow  in  the  north  was  opposing  the 
environment  which  may  be  of  consequence.  Also  with  sufficient 
time  resolution,  the  propagation  of  these  vortices  and  rotors 
may  be  studied.  The  solenoidal  term,  although  believed  to  be 
small,  may  contribute  in  the  head  of  the  gust  front.  This  may 
be  important  in  the  formation  of  the  rotor.  The  elongated 
region  of  horizontal  vorticity  mentioned  for  the  5  August  case 
also  coincides  with  a  perturbation  temperature  gradient. 
Perhaps,  this  is  a  hint  of  possible  solenoidal  activity. 


APPENDIX 


1.  Vorticity  Equations:  By  taking  the  curl  (VX  )  of  the  equa¬ 
tion  of  motion  in  vector  form,  the  following  vorticity  equation 
results : 

dcf  /dt  -  -q  (V-V)  +  q  .  7V  -  7  X  (oVP)  +  7  X  F 
a  a  3 

•* 

where  V  is  the  three  dimensional  wind  vector; 

■<}  »  q  +  225,  is  the  absolute  vorticity  vector; 

a 

d/dt  -  3/3t  +  V  •  7.  +  w3/3z,  is  the  total  derivative; 
h  h 

V  -  7^  +  \d/dz,  is  the  3-D  del  operator; 

q  -  U  +  U  *  ^C,  is  the  relative  vorticity  vector; 

2Q  -  lO  +  j  2Qcos<|>  +  k  2Qsin$,  is  the  angular  velocity 

^  ^  A 

vector  of  the  earth,  i,  j,  and  k  are  unit  vectors  toward  the 
x-,  y-,  and  z-axes,  respectively. 

Equation  (1)  is  the  complete  form  of  the  vorticity  equa¬ 
tion  in  Cartesian  coordinates  in  the  relative  frame  of  refer¬ 
ence.  Where  (1)  is  applied  to  the  convective  scale  motion, 
such  as,  the  meso-(J  scale  (Orlanski,  1975),  the  contribution 
due  to  the  earth's  rotation  to  the  absolute  vorticity,  2Q, 
becomes  sufficiently  small.  Hence  q^  can  be  replaced  by  q.  It 
follows  that  (1)  can  be  rewritten  as: 

dq/dt  *  3q/dt  +  V«Vq  +  w3q/3z 

-  -q(V  •  V. +3w/3z)  +  q«^V  -  V  X  (crVP)  +  V  X  F 
n  n 
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where  term  A  is  the  local  time  rate  of  change  of  the  rela¬ 
tive  vorticity  vector; 

term  B  is  the  horizontal  advection  of  q; 

term  C  is  the  vertical  advection  of  q; 

term  D  is  the  divergence  term; 

term  E  is  the  tilting  term; 

term  F  is  the  solenoid  term; 

and  term  G  is  the  friction  term. 

Term  D  can  be  expanded  as: 

D:  -('lZ  +  jn  +'kC)(3u/3x  +  3v/9y  +  3w/3z)  (3) 

Similarly,  terms  E,  F,  and  G  can  also  be  expanded  as: 

E:  ((i£  +  jn  +  kC)»  V)  (lu  +  jv  +  lew)  (4) 

<5> 

F:  -V  X  (aVP)  -  -  Va  X  VP  ■  -  £( 3a/3y3p/9z-3a/9z3p/3y ) 

+  ^(3a/3z3p/3x-3a/3x3p/3z)  -  IkC 3a/3x3p/3y-3a/3y3p/3x) 

(6) 

G:  (i9/3x  +  ^3/9y  +  1e3/9z)  X  (if*  +  ^Fy  +  icFz>  - 

*i(3F  /3y-3  /3z)  +  3 (3F  /9z-9F  /9x)  +  1c(3F  /9x-3F  /3y) 
z  y  Jx  z  y  x 

Remember,  the  vorticity  equation,  thus  far,  has  been 
operated  in  the  coordinate  system  with  respect  to  the  earth. 
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In  a  convective  storm  with  a  constant  translational  speed,  c, 
it  is  essential  to  use  the  storm-relative  coordinate  system  or 
rather,  moving  coordinates  with  the  system's  origin  at  the 
center  of  the  storm.  This  is  particularly  needed  in  a  thermo¬ 
dynamical  retrieval  study.  Recall,  the  velocity  vector,  V,  is 
the  total  wind  vector,  including  both  the  storm-relative  wind 

a 

vector,  V  ,  and  the  storm  motion,  c,  i.e., 

V  «  V  +  c  (7) 

r 

If  the  storm  motion  is  nearly  stationary,  c**0,  then 

V-V  (8) 

r 

For  the  storm  cases  under  investigation  using  the  JAWS 

data,  c  is  very  small,  therefore,  no  space-time  conversion  is 

needed.  In  other  words,  V  can  be  treated  as  V  ,  etc.  However, 

r 

if  the  system  is  moving,  the  storm-relative  wind  vector, V  , 
must  be  used  to  replace  V.  Remember, 


V-V  -  V 


•  ( V  +  c)  “  V-Vr  + 


V  X  V  -  V  X  (V 


vr  .  C)  - yA  *‘nvr 


(9) 


9u/9x  «  9/9x  (u  +c  )  -  9u  /9x  +  9c//9x  -  9u  /3x 
r  x  r  * 


+  9c/^9: 


and  so  forth. 


From  now  on,  V  can  be  thought  of  as  V^;  however,  for  con¬ 
venience,  the  subscript  r  will  be  dropped.  Furthermore,  the 
friction  term,  in  general,  is  one  order  of  magnitude  smaller 
then  the  other  terms  in  (2)  and  will  be  dropped  hereafter. 
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2.  Vorticity  components  in  moving  coordinates  (x,y,z): 
Three  vorticity  component  equations  in  moving  coordinates 
(x,y,z,t)  can  be  derived  from  (2)  by  performing  dot  (i,j,and  k) 
products  on  (2).  This  yields: 

3£/3t  -  -V,. -w3£/3z  -£(V  •  Vv+3w/3z)  +(£3u/3x+n3u/3y+C3u/3z) 
h  h  h  h 

-(3a/3y3p/3z-3a/3z3p/3y ) 

3n/3t  -  -V  .  V.n  -w3q/3z  -q(V  •  V +3w/3z)  +  (£3v/3x+r|3v/3y+C3v/3z) 
n  n  n  n 

-(3a/3z3p/3x-3a/3x3p/3z)  (10) 

3C/3t  -  -V  •  V^C  -w3C/3z  -C(V  •  Vw  +  3w/dz)  +(£3w/3x+n3w/3y+C9w/3z) 
h  h  n  n 

-(3a/3x3p/3y-3a/3y3p/3x) 


where  £  *  3w/3y-3v/3z  is  the  vorticity  component  about  the  x- 
axis; 

q  *■  3u/3z-3w/3x  is  the  vorticity  component  about  the  y- 
axis;  and 

C  *  3v/3x-3u/3y  is  the  vorticity  component  about  the  z- 

axis . 

3.  Computation  of  vorticity  budgets:  Equation  (10)  can 
be  used  to  compute  te  vorticity  generation  of  a  mesoscale  sys¬ 
tem,  e.g.,  the  mesocyclone,  the  microburst  misocyclone,  etc. 
Excepting  the  solenoidal  term,  which  uses  the  information  from 
the  Doppler  derived  wind  field  and  is  subject  to  larger  uncer¬ 
tainty,  the  other  terms  can  be  determined  from  the  three  wind 
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components  (u,v,w)  and  their  derivatives. 
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